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We study asymptotic dynamics of photons propagating in the polarized vacuum of a locally 
de Sitter Universe. The origin of the vacuum polarization is fluctuations of a massless, 

<N ■ . 

■ minimally coupled, scalar, which we model by the one-loop vacuum polarization tensor of 



scalar electrodynamics. We show that late time dynamics of the electric field on superhorizon 
scales approaches that of an Airy oscillator. The magnetic field amplitude, on the other hand, 
asymptotically approaches a nonvanishing constant (plus an exponentially small oscillatory 
component), which is suppressed with respect to the initial (vacuum) amplitude. This implies 
^ | that the asymptotic photon dynamics is more intricate than that of a massive photon obeying 



the local Proca equation. 
PACS numbers: 98.80. Cq, 04.62. +v 



INTRODUCTION 



Even though we have recently witnessed significant progress in understanding the dynamics of 



gauge fields 



work |ll E we computed the one loop vacuum polarization iron! charged, massless. minimally 
coupled scalars in a locally de Sitter background. We also evaluated the integral of the position- 
space vacuum polarization tensor against a tree order photon wave function and found a nonzero 
result which grows with the same number of scale factors as the mass term of the Proca equation. 



photons) during inflation, our understanding remains quite rudimentary. In earlier 
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What we have so far done has a direct analogue in flat space quantum field theory. Because the 
background possesses spacetime translation invariance it is much simpler to work in momentum 
space. Gauge and Lorentz invariance imply that the vacuum polarization is proportional to the 
transverse projection operator, IF 1 " (p) = (p 2 "rf' v — p^p u )Ii{p 2 ). Here rf' v is the Minkowski metric, 

TV = diag(-l, 1,1,1), tT = diag(-l, 1, 1, 1) , (1) 

and p 2 = p IJ, p u rj f j_ u . One checks that the photon remains massless by evaluating p 2 Ii{jp 2 ) on the 
tree order mass shell, p 2 = 0. A nonzero result — as in the Schwinger model jj] — proves that 
quantum corrections change the tree order mass shell, but it does not determine the fully corrected 
mass shell. For that one must find the pole of the full propagator, that is, solve for p 2 such that, 



P 2 



i-n( P 2 ) =o. (2) 



What we have done for photons during inflation is the analogue of computing II (p 2 ) at one loop 

lp 2 =0 



and showing that p 2 U(p 2 )\ p 2 =0 ^ 0. What we do in this paper is the analogue of solving equation 



Fo introduce the problem, we remind the reader that we work in the conformal coordinate 
system of a locally de Sitter space-time whose metric and inverse are, 

V=«V. <r = a-V- (3) 

Here a is the scale factor, 

«(l) = -^> (r,<0). (4) 

Fhe Fagrangean density of massless minimally coupled scalar electrodynamics in a general metric 
reads, 

Uqed = -\^g w 9 U(T F^F pa - (d M - ieA^*{d v + ieA v )<j> , (5) 

where = d p A u — d^A^. Its reduction to the locally de Sitter background (jSJBJ) is, 

Uqed -> -^7]l»Tf' a F^F pa -a 2 Tr(d»-ieAr)<l>*(d v + ieA v )<l>, (6) 

where we made use of \f—g = a 4 . From this form of the Fagrangean, it is obvious that the photon 
field couples conformally to gravity, with the trivial rescaling, 

Fpv -> F pv (A^ -» A/j) . (7) 
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An important consequence is that tree level photon wave functions take the same form in conformal 
coordinates as they do in flat space, oc e^e lk ' x . 

On the other hand, the minimally coupled scalar field <f> does not couple conformally to gravity, 
which can be easily seen upon the conformal rescaling, <fi — > a _1 in Eq. (jSJ). This has as a 
consequence the extensively studied particle creation during inflation, known as superadiabatic 
amplification. At one loop this conformal symmetry breaking is communicated to the photon 
sector by the cubic and quartic terms in (jfjj). The main objective of this work is to study how this 
changes photon wave functions. 

Provided it couples minimally to gravity, an obvious candidate for the scalar field <ft is the charged 
component of the fundamental Higgs scalar, which at low energies manifests as the longitudinal 
component of the boson. Of course, there may be many more charged scalars of this type 
whose masses lie between that of the Higgs field, mjj ~ 10 2 GeV, and the energy scale of inflation, 
H ~ 10 13 GeV. An important class of such particles would be the supersymmetric scalar partners 
of the standard model quarks and charged leptons. 

In the presence of charged scalar fluctuations the vacuum becomes polarized, such that in 
conformal space times (J3J) Maxwell's equations generalize to Q], 

frrf°dpF<n,(x)+l d 4 x'[^ ct ](x, x')A u (x') = . (8) 

Here [ M n^ ct ](x, x') denotes the retarded vacuum polarization bi-tensor, which has the general form, 

m ct ](x, x') = -[V] Xe (x, x') + [^] 5n\x, x') . (9) 

The two transverse projectors are defined by, 

[^p u \ = [vr^r - rf^Wpdc , rn = ^ ^ v . v - ^ , (10) 

with <9 M = d/dx^, = d/dx'v. The coefficient of [»P U ] in © is, 

Sn 2 (x,x') = Xe(x,x') + * m (x,x'). (11) 

-L + Xm 

where Xe(x,x') and Xm(x,x') denote the electric and magnetic susceptibilities of the vacuum. An 
important consistency check of equations (jSJ)-© is obtained by taking a derivative of Eq. (|SJ). 
The first term vanishes by the antisymmetry of F au , while the integral term vanishes because of 
the transverse structure of the vacuum polarization tensor ©, 

M ["IT] =0 = ^["n"]. (12) 
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This is an immediate consequence of <9 M pP"] = 0, 9„ fP"] = 0, d' u [^P u ] = 0, and S^P"] = 0. 

In a locally de Sitter space-time Q) there is a simple relation between £(x;x'), the invariant 
length from to x /At , and the conformal coordinate interval Ax 2 (x;x'), 



4sin 2 (^He(x;x')J = aa'H 2 Ax 2 (x;x') 



y{x;x 



(13) 



We refer to y(x; x') as the de Sitter length function. The scalar propagator can be expressed in 
terms of y = y{x\ x') and the two scale factors, 

R D-2 



iA(x; x) 



(47t)t 



2 D - 4 r(P>-l)ln(aa')-vrcot(7rf) 



r(f) 



+ E 



1 r(D -l + n)/y 



^ L n T(f+n 



1 + n) /y\n- 



n 



D 



Tin) 



(!)' 



_d n 

2 



(14) 



Note that the homogeneous terms (i.e., y n ) and the constant terms have slightly different propor- 



tionality constants from our previous expression 



HQ. 



This has been done to make ()14|) valid for 



regulating a spaed mw in which D will ultimately be la ken to some dimension oilier than (bur 1^) . 
In D = 4, the propagator (jTl|) reduces to, 



zA(x; x ) 



H 2 jrrrf_ 



iln(P. 2 Ax 2 )-i+ln(2)|. 



(15) 



4vr 2 I Ax 2 2 

In the presence of charged scalar fluctuations described by this propagator, the photon in a 
locally de Sitter space-time sees a polarized vacuum. This effect is characterized by the renormalized 
vacuum polarization bi-tensor, the one loop result for which is 

SB. 

i 



ix rcn 



(x, x') 



a 



32vr 3 



\iipv 



12 



— d 4 ( ln 2 (/j 2 Ax z ) - 21n(^Ax :£ 



+ i 



16tH 



ln(a)£ 4 (x — x') 



1 



+ ln 2 (P/ ,2 Ax 2 ) +21n(P;' 2 Ax 2 

7777' V 



77*77' 



■(ln 2 (Pe 2 



2 /2 , ... 2 Ax 2 )+41n(Pe 2 Ax 2 ) 



(16) 



Here u is the renormalization scale, a = e 2 /47r is the fine structure constant, and H' = 2~ 1 e 1 ^H. 
The polarization tensor (|16|) was obtained using dimensional regularization with the scalar propa- 
gator Q14j). It is fully renormalized, which means that it is an integrable function of x'^ when the 
derivatives are taken after the integration. It is also gauge invariant on account of its transversality 
on both indices (|12|). 

The renormalized vacuum polarization Q16JI is properly an in-out matrix element. To convert 
it into the retarded vacuum polarization [^11^] of Eq. (JHJ) we employ the Schwinger-Keldysh 



formalism as has been described at length elsewhere Q]. The result takes the same form as © 
with electric susceptibility, 

Xe Xelflat Xe| anomaly Xe|de Sitter 

a 



Xelflat = g^9 4 [(ln( M 2 AT 2 )-l)0(A77)0(Ar 2 



967T 2 

Xe I anomaly $ (•£ ) 

a 



XeldeSitter = ^H 2 aa 1 d 2 [(\n(H> 2 At 2 ) + i)0(At/)0(At 2 )] . (17) 

The change in the index of refraction is, 

Sn 2 = ^H A a 2 a' 2 (\ti{H' 2 At 2 ) + 2) 0(At?)0(At 2 ) . (18) 

In these formulae a = a(rf), a' = a(rj'), H' = 2~ 1 e 1 / 4 H, At 2 = (77 — r/') 2 — \\x — x'\\ 2 , A-q = rj — rj', 
and 9{x) = 1 for x > 0, 9{x) = for x < 0. 

The electric susceptibility in (|17|) can be neatly split into three contributions. The first contribu- 
tion Xelflat is conformally invariant. When expressed in terms of conformal time, this contribution 
becomes identical to the one-loop electric susceptibility in Minkowski (flat) space-time. This term 
contains no scale factors, and it is renormalized precisely as in flat space. The second contribution, 
XA~* -ne S tao the trace anotnaly, which arises as a _ce of perfect cancellation 
between the counterterm and the local one-loop diagram in dimensional regularization [2j|. The 
anomalous contribution depends only weakly on the scale factor a. The trace anomaly was first 
considered in the context of electrodynamics in conformal space-times in Ref. [f|, while for a dis- 
cussion of scalar electrodynamics and nonabelian gauge theories see Refs. Q and The third 
contribution to the electric susceptibility, Xelde Sitter m (UH> is purely inflationary, completely finite 
and grows linearly with the scale factors a and a'. And finally, Sn 2 in (|18j) is the second inflationary 
contribution, which grows quadratically with the scale factors a and a' . 

In section [U] we argue that only the one loop terms Xe|de Sitter an d Sn 2 need be kept in equation 
(JHJ). We also give the integro-differential equation which describes a spatial plane wave with arbi- 
trary time dependence. (The derivation is reserved for the Appendix.) In section 0X11 we analyze 
the integro-differential equation for superhorizon photons at asymptotically late times. We show 
analytically that while the vector potential approaches a nonzero constant, its time derivative be- 
haves like an Airy oscillator. Numerical analysis is also presented. In section HVl we translate these 
results into statements about the physical electric and magnetic fields. 
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II. PHOTON DYNAMICS WITH VACUUM POLARIZATION IN INFLATION 

We shall now study the dynamics of the photon field by including only the genuinely new 
inflationary terms contributing at one-loop to the vacuum polarization tensor, Xe|de Sitter 

(|T7|) and 

5n ()18|) . since their contribution to the vacuum polarization grows rapidly with the scale factor. 
Indeed, 5n oc a 2 , and the electric susceptibility, Xelde Sitter °c a. From d a a = 5® Ha 2 it follows 
however, that this term can also give an a 2 in Eq. (jHJ) and has to, therefore, be treated on an 
equal footing as 5n 2 . Moreover, we expect that, when compared with the one-loop contributions 
to 

Xe|de Sitter an d 5n 2 , the contributions of higher loops can be neglected, since they are down by 
more powers of a, but can give no extra powers of a. 

A general solution can be expressed as a linear combination of spatial plane waves in Coulomb 
gauge, 

A fl (x) = e fl (k,7 1 )e i ^, eiki = 0. (19) 

In the Appendix we show that eo(k,r]) vanishes in the absence of sources, just as it does in flat 
space. The Appendix also describes how to perform the integral over x' to obtain the following 
integro-differential equation for the time dependent polarization vector, 

aH 2 



-(ctf + k 2 )ei(k, V ) + 



irk 



(d 2 + k 2 ) Pdrfa^fkAr,, A) ei (k^') 

+ H 2 a 2 dr]'a' 2 siu{kAr})ei{k,r]') =0. (20) 
J no 

Here Yi(x, Q is, 

Yx(x,() = sin(x)[21n(2a;/C) + 2 + ci(2x) - j E - ln(2x)] - cos(x)[si(2x) + vr/2] . (21) 
Upon acting the time derivatives upon the first integral, Eq. (|20[) can be recast as, 
(do + k 2 )ei(k,r]) + a 2 ~~~~~~{ ^ m ( a ) " ~^$,'n) - ^ d V 'a' 2 sin(fcA7/)ei(^ rf) 

a J no Arj J m V rfj 

- ln(l - i) ei (fc,77o)} =0. (22) 

This integral-differential equation is one of our main results. It should describe photon dynamics 
accurately from the beginning of inflation, when 0(770) = 1 and the one loop terms are negligible, 
all the way to late times, when a 3> 1 and the one loop terms become important. Another of 
our main results is that, although certain aspects of this dynamical system can be described in 
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terms of local field equations, the system as a whole is essentially nonlocal. There is solid physics 
behind this: the nonlocal terms reflect the influence of scalar fluctuations within the past light- 
cone. Note that this is perfectly consistent with causality. Note further that, since ln(l — rj/rj') 
diverges logarithmically at 7/ = 77, the last integral cannot be recast in the form of an integral over 
a kernel K = K(t7,t/) multiplied by the polarization vector, H K(r],r]')ei(k,r]'). Instead, its form is 
that of a nonlocal conductivity, J^ o K'(rf,rf)&Qei(k,rf). We shall now address the relevance of these 
nonlocal (potentially dissipative) contributions. 



III. LATE-TIME DYNAMICS OF SUPERHORIZON PHOTONS 



A. Analytical analysis 

We begin by performing a partial integration on Eq. (|22|) . 

{dl + k^(k, r,) + a 2 F 2 V { ( 111 (|) " I + ™) 1) 

vo kr\ - 



TT 



c\(kAr] ) + 



d 'ei(k,T]' 
sm(kAr]o) 



krjo 



(23) 



It is now convenient to rewrite this equation in terms of the number of e-foldings, N = ln(a" 
-H-Hrj), 

(d% + d N + w 2 a- 2 )ei(w, N) + ^ j ( In (w) - ~ + je^Mw, N) 

<1X' ci^fa'" 1 -a- 1 }) 



N 





2a 

TT 



sin (w[a! 1 — a *]) 



wa 



1-1 



d %(w,N' 



ci 



(w(l - 1/a) 



sin 



(w(l - 1/a) 



IV 



€i(w,0) 



(24) 



The new dependent variable is, ei(w,N) = ej(/c,ry), where uv = k/H S> 1. The scale factors are 
a = e N and a' = e N . When N becomes large the terms containing negative powers of a can be 
dropped to give, 



N 



dN' 



ciiwa 



sin (wa 



wa' 



d N >€i(w,N') \ ~ 



2a 



TT 



ci(w) 



sin to 



ei(™,0). (25) 



It is easy to see from equation (|25[) that €i(w,N) approaches a constant plus an exponentially 
small term. First note that the assumption is self consistent. If we make it then, for large enough 
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N', the integrand of the nonlocal term falls off. This means that the integral is independent of its 
upper limit and the constant is determined by the relation, 



In (w)-^ + j^j€i(w,oo)-J^ dN' 



f-u sin [wa! 1 ) n 



Cl TO 



wa 



i-t 



d N ,ei{w,N') 



ci(w) - 



sin(w) 



w 



€i(w,0) . 
(26) 



However, because the integral depends upon d^ie^w, N') for finite N', when it is still significant, 
the actual value of €i(w, oo) must be determined numerically. 

To get the rate at which Ii[w,N) approaches ei[w, oo) we differentiate (|23)l with respect to N. 
Neglecting exponentially small terms gives the following local equation for djsrei[w,N), 



{afa + d N + ^(N-^)}d N ei(w,N)~0. 



(27) 



The general solution can be expressed in terms of Airy functions, 

1 7T 



d N €i(w,N) ~ ci-V^Ca Ai(-(2a/vr) 1 / 3 



N- 



4 8a 



+ C B Bi -(2a/vr) 



1/3 



N- 



1 7T 

4~8a 



}. (28) 



Making use of the asymptotic expansion of the Airy functions Eq. (|28|) can be integrated to 
give asymptotically, 



iV>lm», (29) 



; (C A sin <p\N, a] - C B cos <p\N, a] 



\ ^i/2 ( 2a/vr)Vi2(Ar_i_^) 3/4 



where, 



<p[N,a] 



2 


r2a- 


1/2 


N — 


1 


7T " 


3/2 


7T 


3 


- 7T - 




4 ~ 


8a- 




4 



(30) 



Of course the constants Ca and Cb are fixed by the behavior of €i(w,N) before the asymptotic 
form obtains. They can only be determined numerically. 

B. Numerical analysis 

By performing a numerical study of Eq. 1)22(1 . we now confirm the analytical results of section ITTT1 
above. The following approximation turns out to be quantitatively justified. Upon neglecting the 
latter two integrals, Eq. (|2*2*)) reduces to, 

H 



{d{, + k z )e i (k, V ) + a 



2 2aH 2 rr, . , , , .. 

{[ln(a)- 7 £ + ln( — 



7T 



i(k,r)) - C{k,f])} ~ 0, 



where we define, 



C(fe,r?) = 



h 2 n . , /a 



drj'a' sin(kArj)ei(k,rj') , 



(31) 



(32) 



no 
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N[e-foldings] N[e-foldings] 

FIG. 1: A comparison of the exact (numerical) solution of the integro-differential equation (1221) with that 
of Eqs. jUBEnil for k = 40H. At late times the photon amplitude, ej, approaches a constant. The difference 
between the two solutions is smaller or of the order the thickness of the lines, and quantitatively, it reaches 
about 0.7%. 

The function £(k,<r]) obeys the following differential equation and initial conditions, 

(^2 +k 2 )C = a 2 H 2 ei C(Mo) = 0, — C(fc, r,)^ = . (33) 

To study how quantitative this approximation is, in figure ^ we compare the exact (numerical) 
solution of Eq. (|22j) with the approximate one obtained by solving Eqs. (|31H33|) . As can be seen 
in figure ^ the agreement is quite good. This implies that the two final integrals in (|22[) are 
dynamically irrelevant. 

Figure n shows that, at asymptotically late times, €i(w,N) is dominated by the constant value 
€i(w, oo), whose real part is slightly suppressed with respect to its classical asymptotic value of 
cos (40) ~ —.67. In figure [21 we subtract the asymptotic constant from the solution of Eqs. (|31H33|) . 
multiply by e N > 2 , and compare the result with the Airy function. The figure displays an oscillatory, 
slightly damped profile, in reasonably good agreement with our result Q29j) . 

IV. DISCUSSION 

It is illuminating to translate our results for the vector potential, Ai(x) = €j(fe, r])e tk ' x , into 
statements about the behavior of the electric and magnetic fields. In Coulomb gauge and using 
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FIG. 2: A comparison of the solution of H31H33|I with that of the Airy oscillator 1)27(1 . 

conformal coordinates, the densities of electric and magnetic field lines per physical 3- volume are 
given by the following expressions, 



E\x) 



d Ai(x) 



B\x) 



e ijk djA k (a 



a 2 (rf) a 2 (r]) 
For spatial plane waves with o(ry) = e N and w = k/H we have, 



(34) 



'E\x) = He- N d N ei{w,N) 



-ik-x 



B\x) 



-it 



ijk k j e- m e k (w,N) . (35) 



Keeping track only of powers of the scale factor, we can summarize the results of section II I II by 
saying that the full quantum system results in the following asymptotic behavior, 



rrrTd P F lTV (x)+ dV^lM*. x')A„(x') = 



E l ~ a 2 and B % ~ a 



(36) 



It is interesting to contrast the asymptotic behavior (|36|) of the full, nonlocal system with that 
of various local models. The tree order system gives, 



rrrTdpF^O 



E i ~ a~ 2 and B i ~ a' 2 



(37) 



Although the behavior of the magnetic fields agree as regards powers of a, we saw in Figure [J that 
the magnetic field of the actual system analyzed in section III B is smaller by a factor of about 
2/3. On the other hand, the electric field of the actual system is enhanced by an enormous factor 

of 0,2. 

Adding a Proca term with fixed mass m 7 = gH results in the following asymptotic behavior, 



(r] pa d p F a 



g 2 H 2 a 2 A u )=0 



E i ~ a-i+sV 1 "^ 2 and B l ~ aT* + S 



(38) 
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The best perturbative fit to the actual system has g 2 = 2aln(u;)/7r P, [J but we see that the 
agreement between (j3fi)) and (|38jl is not good. The constant mass system has too much electric 
field and not enough magnetic field. 

If the mass term grows as m 2 = g 2 H 2 In (a) one finds, 



V 



7 

jiv I pa 



(ri pa d p F au -g 2 H 2 ln(a)a 2 A^=0 => E i ~ a~l and B* ~ o~a . (39) 



This model is suggested by mean field theory 
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11] with g 2 = 2a /it. It gives good agreement 



for the electric field but results in too little magnetic field. In particular note that, whereas the 
actual polarization vector approaches a nonzero constant, e» (10,00), the polarization vector of (|39|) 
approaches zero. 

While none of the local models described above gives a perfect fit, the effect of vacuum po- 
larization is to vastly enhance the asymptotic electric field with respect to its classical value in a 
manner most nearly described by the growing mass of (|39j) . Without regard to local analogues, it is 
very apparent that vacuum polarization alters the kinematical properties of superhorizon photons 
during inflation. On the other hand, the fact that the asymptotic polarization vector is some- 
what suppressed from its classical value seems to indicate that there is no significant production 
of photons. 

Exactly the opposite results seem to pertain for massless fermions which are Yukawa coupled to 
a massless, minimally coupled scalar [121 ] . In that case chirality conservation protects the particle 
from gaining a mass, but the amplitude of the wave function grows faster than exponentially during 
inflation. The physical explanation seems to be copious particle production. 
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Appendix: Reduction of the polarization integral 

Here we show some details of the calculation of the inflationary contribution to the vacuum 
polarization integral in (jHJ) with the general plane wave (|19|) for the photon field, A v (x') = 
e u {Ki)^*' (hei = = e ). 

We begin our analysis by noting a useful identity, 

- p*p"] oaf = - aa'^o V 2 + rf >a'(V 2 - $}) + H 2 a 2 a' 2 (l - A V d )] 

- Haa ,2 5° p d u + Ha 2 a'3»5 u + aa'(<9^ + 6%8 u )d - aa'B^d u . (40) 

In deriving this we made use of the relations, 

d' p aa' = aa'd' p + Haa' 2 5° p and d a aa' = aa'd a + Ha 2 a'5° , (41) 

which imply, 

d' p d a aa' = H 2 a 2 a' 2 5° p 5° a + Haa' 2 5 p d a - Ha 2 a'd' p 5° a - aa'd' p d a . (42) 
Note also that we employ the notation, 

3» = rf % di , t v = *f v + ^o^o • (43) 

When a derivative operator stands on the right of aa! , then it acts on a function of x M — x >fi and we 
can write d p = —d p . Moreover, due to the (spatial) transversality of the photon field, d v A u {x) = 0, 
the terms in (|4Ujl containing d v vanish identically. 

Consider first the 0-th component of the photon field equation (jSj). 

a'A)(x') = 0, (44) 

where we made use of Eqs. (|17I18|) and l|40|l. and of the Coulomb gauge condition, V • A = 0. Note 
that only the first term in (|40|) contributes to (|44|). Furthermore, due to the spatial transverse 
structure of the Lorentz breaking term, [^P^Sn 2 , it does not contribute to ()44|) . 

At 7] = r] the integral drops out of and we conclude that V 2 A (x, r] ) = 0. We can make 
use of residual gauge freedom to make the unique solution be Aq(x, t]q) = 0. The same thing can be 
done for the first conformal time derivative. But then the equation implies that Aq(x, r)) vanishes 
for all r\. Therefore, it is consistent to assume eo(k,r]) = 0. This discussion also implies that the 
terms in (j4U|) which contain 5® cannot contribute to the photon dynamical equation (JHJ). 



V 2 A n 



an 



V 2 d 2 / (fx' ' 8(^)6 '(At 2 ) ln(F /z Ar 2 ) + 1 
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Consider now the spatial components of equation (JHJ). We start with the de Sitter contribution to 
the electric susceptibility, Xe|de Sitter <|17|L which with the help of (|4T)|) (only the terms proportional 
to rf 11 contribute), 



mtl /cfx'[^]x e (x,x')|de Sitter A V { X ') 



" R )ae rk S (d 2 + k 2 ) 2 f drj'a^kAr,, ^forf) 



2nk 3 



'/o 



gff 4 
2irk 3< 



-V fc -(9 2 + k 2 ) I d V 'a' z (l - Ar 1 d )E 1 (kA V , -^j)e t (k, V ') . (45) 



/ /2 



'/ii 



k 



H' 



Here r/Q = —1/H denotes the conformal time at the beginning of inflation, at which a(ryo) = 1- We 
define k = \\k\\, Ar] = r/ — rf. We also made use of the following elementary integrals, 

rAr) 

Jo 

= 4tt / A Vdr^Mfln \h' 2 (A V 2 - \\Ax\\ 2 )} + l\ 
Jo kr I- L J J 



d 3 Axe- ilA *{ In [tf'^A?? 2 - ||Ax|| 2 )] + l} 



—Ar] 2 J xdxsin (kAr]x){ ln(l - x 2 ) + 2ln(H'Arj) + l} 



47T 



Zi(kAri,k/H'), 



(46) 



where Ax = x — x' 



ciu 



H n (z,C) = [sin(z) -xcos(z)][21n(;z/C) +n] + (*) 

z 2 £(z) = z 2 [ z'dz'sm(zz')\n(l-z' 2 ) 
Jo 

= 2sin(z) - [cos(;z) + zsin(z)][si(2z) + it/ 2] + [sin(z) - z cos(z)][ci(2z) - j E - ln(z/2)] 
f°° sinfi) f z sin(i) , 7r 



cos(t) 



dt 

t jo 



cos(t) — 1 



and <y E = -^(1) = -{d/ dz)\u[Y{z) 



t 



\z=l 



dt + ~f E + ln(z) 



(47) 



0.577 215 664.. is Euler's constant. Now making use of, 



Yi(z,() = -(^ + l)3i(*,0 



(48) 



= sin(z) [2 ln(2z/C) + 2 + ci(2z) - j E - \a(2z)] - cos(z) [si(2z) + tt/2] 
(«a z -l)(Sg + l)3i(ar,C) = 4sin(z)-2H 2 (z,C), (49) 



the integral (|43|) reduces to 

-Vi„J d 4 x'^P u ]xe(x,x')\ dcSittev A u (x') 

aH 4 



+ ■ 



irk 



—ae ik - x (d 2 + k 2 ) J drfo/Y^A^—^Krf) 
a 2 e its [ V drj'a ,2 \2sm(kAr])-E 2 (kAr), e^k, rf) . (50) 
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The contribution from 5n 2 to the spatial components of Eq. © can be obtained quite easily by 
the methods analogous to those employed for the contribution from Xe I de Sitter; the result can be 
written as, 

■a~e"~ " 

'Vo 



mil d 4 x'[^]5n 2 {x : x>) A v {x') 



irk 



_ a 2 e ik.s f V dr j' a ' 2 E 2 (kA v , AW£,t/). (51) 

Jr?n v ti I 



Remarkably, this contribution is fully canceled by the analogous one in (J5(J|). Keeping the terms 
that contribute at order a 2 , we have, 

mil f d A x'[^ et ]{x,x') A v (x') = ^ / d 4 x'(-[^]xe(x,x')\ deSiUcv + [^}5n\x,x'))A u (x')+0(a) 



(52) 



" H2 . a ^{dl + k 2 ) f V dri'a'Y, (kAr,, A) ei (k, rf) 



wk 
2aH 4 



irk 



+ —a 2 e lk - x / dr]'a' z sm(kArj)ei(k, rf) + 0{a) , 



'/o 



where 0(a) indicates that the neglected terms can contribute to the vacuum polarization at most 
linearly in a. 

Now from the small argument expansions of the sine and cosine integral (c/. Eq. [8.232] in 
Ref. Q), 

-l) n+1 z 2n_1 



si(z) 



7T 



+ E 



Cl 



2 ^(2n-l)(2n-l)!' 

-l) n z 2n 



(z) = 7 £ + ln(z) + ^ 



£i(2n) (2n)!' 

we obtain easily the following small argument expansion of Y\ , 



\z\ < 1 
M < 1 



(53) 



Y 1 (z,C) = 2zln(^)-Iz 3 [ln(^)-^]+0(. 5 ln(.),z 5 ). 



C 



2z\ 31 



(54) 



C / 12- 

This implies that, attempting to act with the derivative d$ on the integral in 1)52(1 . would result 
in a divergent contribution at the upper limit of integration, rf = rj. In order to overcome this 
difficulty, observe first that one time derivative may be taken, 



^ jy v 'a'Y 1 (kAr l ,^)e i (k,r l ') = J* drj'a' [d z Y x (z, A) 



z=kAr) 



ei(k, rf) . 



(55) 



Note we can write, 



^Yi(z,C) = cos(z)[21n(2z/C) + 2 + ci(2z) - j E - ln(2z)] + sin(z)[si(2z) + vr/2] 

= d z Y x (z, C)-2 ln(2z/C) -2] +2 ln(2z/C) + 2 , (56) 
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where in the last line we have added and subtracted the troublesome term which is logarithmically 
divergent at the upper limit of integration, r( — * r). We can now act with the second time derivative 
on the first term, to arrive at, 

d o +k2 r^ v ji.A„ k \,jt^ _ n r^/ 1 - 008 ^). 



k 



v 



At? 



+ 2d drj'a' ln(2F'Ar?) + 1 e;(fc,r/) 



(57) 



Now from, 



a' ln(2i?'Ary) 



1 



#77' 



H-*«V)-£i$* 

n=l ' 



imV2W) + i;A(^' 

2 ^— ^ n z \77 ' 

n=l 1 



we can rewrite the second integral in (|57jl as follows, 



dr/'a' ln(2^ / Ar ? )e i (fc, 77') 



1 



5ln 2 (-2ilV) + E^(5)1^(^V) 



n=l 
oo 



- / ^/> 2 (-2W) + Ei(5) B 



— (1 

71=1 ' 



77'= J? 

v'=vo 



dr) 



7 ei(k,r/). (58) 



Its time derivative equals, 



d / A/a' ln(2i/'A?7) + 1 ei(k,rf) 



+ a /'Win (l-^)^(fc, V) 



ln(-2F'? ? ) + 1 

r77 ^ 

dr/' In ( 1 - 
'no V »T y a»? 

Upon collecting all of the terms, Eq. ((521) reduces to, 

,2atf 2 



(£,77) + a In (l - — Wfc,77o) 



(59) 



T^riV^t]^')^') 



7T 



# 2 /"'? 

~k 



ik-x 



drj'a' 2 sm(kArj)ei(k, rf) 



iiit 



1 n 1 , .1- cosfkArj) , r ., 
a At? 

-/>' b ( 1 -^^<A"'» 

- In (l - 770) , 



(60) 



where we took account of 770 = —1/H. When o > 1, the last local term, which depends on the 
initial photon amplitude ei(k,r/o), contributes as 0(a~ 1 ), such that, in the limit when a 1, it 
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can be consistently neglected. 
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